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ABSTRACT: 

In this paper.we describe some central mathematical problems in medical imaging. The subject has been undergoing 
rapid changes driven by better hardware and software. Much of the software is based on novel methods utilizing 
geometric partial differential equations in conjunction with standard signal/image processing techniques as well as 
computer graphics facilitating man/machine interactions. As part of this enterprise. researchers have been trying to 
base biomedical engineering principles on rigorous mathematical foundations for the development of software 
methods to he integrated into complete therapy delivery systems. These systems support lite more effective delivery of 
many image-guided procedures such as radiation therapy. biopsy. and minimally invasive surgery. We will show how 
mathematics may impact some of the main problems in this area including image enhancement, registration. and 
segmentation. This paper [1] describes image processing techniques for Diffusion Tensor Magnetic Resonance. In 
Diffusion Tensor MRI,a tensor describing local water diffusion is acquired for each voxel. The geometric nature of 
the diffusion tensors can quantitatively characterize the local structure in tissues such as bone.inuscles,and white 
matter of the bra.in.The close relationship between local image structure and apparent diffusion makes this image 
modality very interesting for medical image analysis. We present a decomposition of the diffusion tensor based on its 
symmetry properties resulting in useful measures describing the geometry of the diffusion ellipsoid. A simple 
anisotropy measure follows naturally from this analysis. We describe how the geometry. or shape. of the tensor can be 
visualized using a coloring scheme based on the derived shape measures. We show how filtering of the tensor data of 
a human brain can provide a description of macroslructiiral diffusion which can be used for measures of fiber-tract 
organization. We also describe how tracking of white matter tracts can be implemented using the introduced 
methods. 
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I. INTRODUCTION 

Medical imaging has been undergoing a revolution in the past decade with the advent of faster, more 
accurate, and less invasive devices. This has driven the need for corresponding software development which in 
turn has provided a major impetus for new algorithms in signal and image processing. Many of these algorithms 
are based on partial differential equations and curvature driven flows which will be the main topics of this 
survey paper.Mathematical models are the foundation of biomedical computing. Basing those models on data 
extracted from images continues to be a fundamental technique for achieving scientific progress in 
experimental,clinical,biomedical,and behavioral research. Today,medical images are acquired by a range of 
techniques across all biological scales, which go far beyond the visible light photographs and microscope images 
of the early 20th century. Modern medical images may be considered to be geometrically arranged arrays of data 
samples which quantify such diverse physical phenomena as the time variation of hemoglobin deoxygenation 
during neuronal metabolism,or the diffusion of water molecules through and within tissue. The broadening scope 
of imaging as a way to organize our observations of the biophysical world has led to a dramatic increase in our 
ability to apply new processing techniques and to combine multiple channels of data into sophisticated and 
complex mathematical models of physiological function and dysfunction. A key research area is the formulation 
of biomedical engineering principles based on rigorous mathematical foundations in order to develop general - 
purpose software methods that can be integrated into complete therapy delivery systems. Such systems support 
the more effective delivery of many image-guided procedures such as biopsy,minimally invasive surgery,and 
radiation therapy. In order to understand the extensive role of imaging in the therapeutic process, and to 
appreciate the current usage of images before,during,and after treatment, we focus our analysis on four main 
components of image-guided therapy (IGT) and image guided surgery (IGS): localization, targeting, monitoring, 
and control. Specifically, in medical imaging we have four key problems: 
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(l)Segmentation- automated methods that create patient-specific models of relevant anatomy from images;(2) 
Registration - automated methods that align multiple data sets with each other;(3) Visualization - the 
technological environment in which image-guided procedures can be displayed;(4) Simulation - softwares that 
can be used to rehearse and plan procedures,evaluate access strategies, and simulate planned treatments. In this 
paper, we will only consider the first two problem areas. However,it is essential to note that in modern medical 
imaging, we need to integrate these technologies into complete and coherent image guided therapy delivery 
systems and validate these integrated systems using performance measures established in particular application 
areas. We should note that in this survey we touch only upon those aspects of the mathematics of medical 
imaging reflecting the personal tastes (and prejudices) of the authors. Indeed, we do not discuss a number of very 
important techniques such as wavelets, which have had a significant impact on imaging and signal processing; 
see [60] and the references therein. Several articles and books are available which describe various 
mathematical aspects of imaging processing such as [67] (segmentation), [83] (curve evolution), and [71, 87] 
(level set methods). Finally,it is extremely important to note that all the mathematical algorithms which we 
sketch lead to interactive procedures. This means that in each case there is a human user in the loop (typically a 
clinical radiologist) who is the ultimate judge of the utility of the procedure, and who tunes the parameters either 
on or off-line. Nevertheless, there is a major need for further mathematical techniques which lead to more 
automatic and easier to use medical procedures. We hope that this paper may facilitate a dialogue between the 
mathematical and medical imaging communities. Diffusion Tensor Magnetic Resonance Imaging (DT-MRI) is a 
recent MR imaging modality. In Diffusion Tensor MRI, a tensor describing local water diffusion is acquired for 
each voxel. Diffusion in tissue can be anisotropic depending on the characteristics of the tissue. For example in 
the white matter fiber tracts the diffusion is mainly in the direction of the fibers. In areas with fluid, such in the 
CSF filled ventricles, the diffusion is spherical, i.e. isotropic. The advent of robust diffusion tensor imaging 
techniques has prompted the development of quantitative measures for describing the diffusion anisotropy. A 
good review by Basser and Pierpaoli can be found in [1]. Since MRI methods in general always obtain a 
macroscopic measure of a microscopic quantity which necessarily entails intravoxel averaging, the voxel 
dimensions influence the measured diffusion tensor at any particular location in the brain. Factors which would 
affect the shape of the apparent diffusion tensor (shape of the diffusion ellipsoid) in the white matter include the 
density of fibers, the degree of myelination, the average fiber diameter and the directional similarity of the fibers 
in the voxel. The geometric nature of the measured diffusion tensor within a voxel is thus a meaningful measure 
of fiber tract organization. With current conventional proton magnetic resonance imaging (MRI) techniques, the 
white matter of the brain appears to be a remarkably homogeneous tissue without any suggestion of the complex 
arrangement of fiber tracts. Although the individual axons and the surrounding myelin sheaths cannot be 
revealed with the limited spatial resolution of in vivo imaging, distinct bands of white matter fibers with parallel 
orientation may be distinguished from others running in different directions if MRI techniques are sensitized to 
water diffusion and the preferred direction of diffusion is determined. Water diffusion in tissue due to Brownian 
motion is random but some structural characteristics of tissues may limit diffusion.In the white matter, the 
mobility of the water is restricted in the directions perpendicular to the axons which are oriented along the fiber 
tracts. This anisotropic diffusion is due to the presence of tightly packed multiple myelin membranes 
encompassing the axon. Myelination is not essential for diffusion anisotropy of nerves as shown in studies of 
nonmyelinated garfish olfactory nerves [3] and anisotropy exists in brains of neonates before the histological 
appearance of myelin [16] but myelin is widely assumed to be the major barrier to diffusion in myelinated fiber 
tracts. Therefore the demonstration of anisotropic diffusion in brain by magnetic resonance has opened the way 
to explore noninvasively the structural anatomy of the white matter in vivo [8, 4, 1, 10]. 



2.1. Generalities 

In 1895, Roentgen discovered X-rays and pioneered medical imaging. His initial publication [82] 
contained a radiograph (i.e. an X-ray generated photograph) of Mrs. Roentgen's hand; see Figure 1(a). For the 
first time, it was possible to visualize non-invasively (i.e., not through surgery) the interior of the human 
body.The discovery was widely publicized in the popular press and an "Xray mania" immediately seized Europe 
and the United States [30, 47]. Within only a few months, public demonstrations were organized, commercial 
ventures created and innumerable medical applications investigated; The field of radiography was born with a 
bang!Today,medical imaging is a routine and essential part of medicine. Pathologies can be observed directly 
rather than inferred from symptoms. For example, a physician can non-invasively monitor the healing of 
damaged tissue or the growth of a brain tumor, and determine an appropriate medical response. Medical 
imaging techniques can also be used when planning or even while performing surgery. For 
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(a) First radiograph of Mrs. Roentgen's hand. 
Fig:l. X-ray radiography at the end of the 19th century. 

example,a neurosurgeon can determine the "best" path in which to insert a needle,and then verify in real time its 
position as it is being inserted. 

III. MAGNETIC RESONANCE IMAGING 

This technique relies on the relaxation properties of magnetically -excited hydrogen nuclei of water 
molecules in the body. The patient under study is briefly exposed to a burst of radio -frequency energy, which,in 
the presence of a magnetic field,puts the nuclei in an elevated energy state. As the molecules undergo their 
normal, microscopic tumbling, they shed this energy into their surroundings, in a process referred to as 
relaxation. Images are created from the difference in relaxation rates in different tissues. This technique was 
initially known as nuclear magnetic resonance (NMR) but the term "nuclear" was removed to avoid any 
association with nuclear radiation.MRI utilizes strong magnetic fields and non-ionizing radiation in the radio 
frequency range, and according to current medical knowledge,is harmless to patients. Another advantage of MRI 
is that soft tissue contrast is much better than with X-rays leading to higher-quality images, especially in brain 
and spinal cord scans. See Figure 2(a). Refinements have been developed such as functional MRI (fMRI) that 
measures temporal variations (e.g., for detection of neural activity), and diffusion MRI that measures the 
diffusion of water molecules in anisotropic tissues such as white matter in the brain. 




Fig: 2(a). Magnetic Resonance Image(brain, 2D axial slice). 

IV. POSITRON EMISSION TOMOGRAPHY 

The patient is injected with radioactive isotopes that emit particles called positrons (anti- 
electrons). When a positron meets an electron,the collision produces a pair of gamma ray photons having the 
same energy but moving in opposite directions. From the position and delay between the photon pair on a 
receptor, the origin of the photons can be determined. While MRI and CT can only detect anatomical 
changes, PET is a functional modality that can be used to visualize pathologies at the much finer molecular 
level.Ths is achieved by employing radioisotopes that have different rates of intake for different tissues. For 
example,the change of regional blood flow in various anatomical structures (as a measure of the injected 
positron emitter) can be visualized and relatively quantified. Since the patient has to be injected with radioactive 
material,PET is relatively invasive. The radiation dose however is similar to a CT scan. Image resolution may be 
poor and major preprocessing may be necessary.See Figure 3(a). 




Figure 3(a). Positron Emission Tomography(brain, 2D axial slice). 
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V. ALGORITHMS & PDES 

Many mathematical approaches have been investigated for applications in artificial vision(e.g. fractals 
and self-similarity,wavelets,pattern theory, stochastic point process,random graph theory; see [42]). In 
particular,methods based on partial differential equations (PDEs) have been extremely popular in the past few 
years [20, 35]. Here we briefly outline the major concepts involved in using PDEs for image processing. As 
explained in detail in [17], one can think of an image as a map I : D — > C, i.e., to any point x in the domain D, I 
associates a "color" I(x) in a color space C. For ease of presentation we will mainly restrict ourselves to the case 
of a twodimensional gray scale image which we can think of as a function from a domain D = [0, 1] x [0, 1] to 
the unit interval C = [0, l].The algorithms all involve solving the initial value problem for some PDE for a given 
amount of time. The solution to this PDE can be either the image itself at different stages of modification, or 
some other object (such as a closed curve delineating object boundaries) whose evolution is driven by the 
image. For example, introducing an artificial time t, the image can be deformed according to 



"* 0) 

where I(x, t) : D x [0, T ) — > C is the evolving image,F is an operator which characterizes the given algorithm, 
and the initial condition is the input image 10. The processed image is the solution I(x, t) of the differential 
equation at time t.The operator F usually is a differential operator, although its dependence on I may also be 

nonlocal. Similarly,one can evolve a closed curve representing the boundaries of some planar shape ( ^ need not 
be connected and could have several components). In this case, the operator F specifies the normal velocity of 

the curve that it deforms. In many cases this normal velocity is a function of the curvature k of ^,and of the 

image I evaluated on ^ .A flow of the form 



is obtained, where N is the unit normal to the curve ^ .Very often,the deformation is obtained as the steepest 
descent for some energy functional. For example,the energy 



(3) 

and its associated steepest descent, the heat equation, 



di_ 



...(4) 



correspond to the classical Gaussian smoothing. The use of PDEs allows for the modelling of the crucial but 
poorly understood interactions between top-down and bottom-up vision 5. In a variational framework, for 
example, an energy £ is defined globally while the corresponding operator F will influence the image 
locally. Algorithms defined in terms of PDEs treat images as continuous rather than discrete objects. This 
simplifies the formalism, which becomes grid independent. On the other hand models based on nonlinear PDEs 
may be much harder to analyze and implement rigorously. 

VI. IMAGING PROBLEMS 

Medical images typically suffer from one or more of the following imperfections: 

• low resolution (in the spatial and spectral domains); 

• high level of noise; 

• low contrast; 

• geometric deformations; 

• presence of imaging artifacts. 
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These imperfections can be inherent to the imaging modality (e.g., X-rays offer low contrast for soft 
tissues,ultrasound produces very noisy images,and metallic implants will cause imaging artifacts in MRI) or the result of a 
deliberate trade-off during acquisition.For example,finer spatial sampling may be obtained through a longer acquisition 
time.However that would also increase the probability of patient movement and thus blurring.In this paper, we will only be 
interested in the processing and analysis of images and we will not be concerned with the challenging problem of designing 
optimal procedures for their acquisition.Several tasks can be performed(semi)-automatically to support the eye-brain system 
of medical practitioners. Smoothing is the problem of simplifying the image while retaining important information. 
Registration is the problem of fusing images of the same region acquired from different modalities or putting in 
correspondence images of one patient at different times or of different patients. Finally,segmentation is the problem of 
isolating anatomical structures for quantitative shape analysis or visualization.The ideal clinical application should be 
fast,robust with regards to image imperfections, simple to use,and as automatic as possible. The ultimate goal of artificial 
vision is to imitate human vision,which is intrinsically subjective.The technique we present below are applied to two- 
dimensional grayscale images. The majority of fhem,however,can be extended to higher dimensions(e.g., vector-valued 
volumetric images). 

6.1. Image Segmentation 

When looking at an image, a human observer cannot help seeing structures which often may be identified with objects. 
However,digital images as the raw retinal input of local intensities are not structured. Segmentation is the process of creating 
a structured visual representation from an unstructured one.The problem was first studied in the 1920's by psychologists of 
the Gestalt school (see Kohler [54] and the references therein),and later by psychophysicists [49, 95]. In its modern 
formulation,image segmentation is the problem of partitioning an image into homogeneous regions that are semantically 
meaningful,i.e.,fhat correspond to objects we can identify.Segmentation is not concerned with actually determining what the 
partitions are. In that sense,it is a lower level problem than object recognition.In the context of medical imaging,fhese 
regions have to be anatomically meaningful.A typical example is partitioning a MRI image of the brain into the white and 
gray matter. Since it replaces continuous intensities with discrete labels,segmentation can be seen as an extreme form of 
smoothing/information reduction. Segmentation is also related to registration in the sense that if an atlas can be perfectly 
registered to a dataset at hand,then the registered atlas labels are the segmentation. Segmentation is useful for visualization,it 
allows for quantitative shape analysis,and provides an indispensable anatomical framework for virtually any subsequent 
automatic analysis. Indeed, segmentation is perhaps the central problem of artificial vision,and accordingly many approaches 
have been proposed (for a nice survey of modern segmentation methods,see the monograph [67]).There are basically two 
dual approaches. In the first, one can start by considering the whole image to be the object of interest,and then refine this 
initial guess. These "split and merge"techniques can be thought of as somewhat analogous to the top-down processes of 
human vision. In the other approach,one starts from one point assumed to be inside the object,and adds other points until the 
region encompasses the object.Those are the "region growing"techniques and bear some resemblance to the bottom-up 
processes of biological vision.The dual problem to segmentation is that of determining the boundaries of the segmented 
homogeneous regions. This approach has been popular for some time since it allows one to build upon the well-investigated 
problem of edge detection(Section 6.2).Difficulties arise with this approach because noise can be responsible for spurious 
edges. Another major difficulty is that local edges need to be connected into topologically correct region boundaries.To 
address these issues,it was proposed to set the topology of the boundary to that of a sphere and then deform the geometry in a 
variational framework to match the edges.In 2D,the boundary is a closed curve and this approach was named snakes. 
Improvements of the technique include geometric active contours and conformal active contours.All these techniques are 
generically referred to as active contours.Finally,as described in [67], most segmentation methods can be set in the elegant 
mathematical framework proposed by Mumford and Shah [69] . 

6.2. Edge Detectors 

Consider the ideal case of a bright object on a dark background.The physical object is represented by its projections on 
the image I.The characteristic function 1 of the object is the ideal segmentation, and since the object is contrasted on the 
background, the variations of the intensity I are large on the boundary flO .It is therefore natural to characterize the 
boundary d as the locus of points where the norm of the gradient |I| is large.In fact, if d is piecewise smooth then III is a 
singular measure whose support is exactly flO .This is the approach taken in the 60's and 70's by Roberts [81] and Sobel 
[91] who proposed slightly different discrete convolution masks to approximate the gradient of digital images. Disadvantages 
with these approaches are that edges are not precisely localized, and may be corrupted by noise. See Figure 4(b) is the result 
of a Sobel edge detector on a medical image.Note the thickness of the boundary of the heart ventricle as well as the presence 
of "spurious edges" due to noise. Canny [14] proposed to add a smoothing pre-processing step (to reduce the influence of the 
noise) as well as a thinning post-processing phase (to ensure that the edges are uniquely localized). See [26] for a survey and 
evaluation of edge detectors using gradient techniques. A slightly different approach initially motivated by psychophysics 

was proposed by Marr and Hildreth [62, 61] where edges are defined as the zeros of + 1 , the Laplacian of a smooth 

version of the image. One can give a heuristic justification by assuming that the edges are smooth curves; more 
precisely,assume that near 
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(a) Original image b) Sobel edge detector 



c) Marr edge detector. 



Fig. 4 Result of two edge detectors on a heart MRI image. 



an edge the image is of the form 



r(x) = ^) 



.(5) 



where S is a smooth function which vanishes on the edge, e is a small parameter proportional to the width of the 
edge, and R — ► [0, 1] is a smooth increasing function with limits. 



In this work we applied a modified version of the recently proposed Line Scan Diffusion Imaging 
(LSDI) technique. This method, like the commonly used diffusionsensitized,ultrafast, echo-planar imaging (EPI) 
technique [12] is relatively insensitive to bulk motion and physiologic pulsations of vascular origin. But unlike 
EPI,LSDI exhibits minimal image distortion,does not require cardiac gating,head restraints or post-processing 
image correction, and can be implemented without specialized hardware on all standard MRI scanners. Here, we 
present a quantitative characterization of the geometric nature of the diffusiontensors,a method for 
characterization of marcostructural diffusion properties, and a display method for showing clear and detailed in 
vivo images of human white mattertracts.The orientation and distribution of most of the known major fiber 
tracts can be identified using these methods. 

7.1. Imaging Parameters 

Suppose, some data were acquired at a hospital on a GE Signa 1.5 Tesla,Horizon Echospeed 5.6 system 
with standard 2.2 Gauss/cm field gradients. The time required for acquisition of the diffusion tensor data for one 
slice was 1 min;no averaging was performed. Imaging parameters were:effective TR=2.4 s, TE=65 ' J i-i>- 1: = 750 
s/bi„ w =5 s/mm^field of view 22 cm, effective voxel size 4.8_1.6_1.5 mm3,6 kHz readout bandwidth, 
acquisition matrix 128_128.The gradient cycle in the LSDI interleaving scheme was modified to provide 
acquisition of more gradient directions and to allow elimination of the crusher gradients. Instead of alternating 
merely between high and low gradient strengths, the modified sequence cycled through eight configurations of 
the diffusion gradients. In all other respects it was identical to the sequence described in [7]. 

7.2. Calculation of Tensors 

For each slice, eight images are collected with different diffusion weightings and noncollinear gradient 
directions. If S I represents the signal intensity in the absence of a diffusion-sensitizing field gradient and V the 
signal intensity in the presence of gradient ^ — W x ' & f> 9z } jme equation for the loss in signal intensity due to 
diffusion is given by the Stejskal -Tanner formula: 




Fig. 4.1. 
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MATERIALS AND METHODS 
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ln(S) = h($>) "'rJ f5/3)g T Dg, 



where " i is the gyromagnetic ratio of hydrogen 1H (protons), L is the duration of the diffusion sensitizing 
gradient pulses and — ^is the time between the centers of the two gradient pulses. The eight images provide 
eight equations for S in each voxel which are solved in a least-squares sense for the 6+1 unknowns: the six 
independent components of the symmetric diffusion tensor, D ,and h i .In the LSDI sequence, it is easy to show 
that cross terms between the slice select gradient for the 180 pulse and the diffusion sensitizing gradients 
account for less than 0.1% of the diffusion weighting,and have therefore been neglected here. Diffusion 

attenuation due to imaging gradients is already factored into ,as is T2 weighting. 

7.3. Geometrical Measures of Diffusion 

In order to relate the measure of diffusion anisotropy to the structural geometry of the tissue a mathematical 
description of diffusion tensors and their quantification is necessary. First, a complete diffusion tensor, D, is 
calculated for each voxel. Using the symmetry properties of the diffusion ellipsoid we decomposed the diffusion 
tensor, and from the tensor basis assigned scalar measures, describing the linearity and the anisotropy, to each 
voxel. The diffusion tensor can be visualized using an ellipsoid where the principal axes correspond to the 

directions of the eigenvector system. Let ^ 1 — ^- — ^ a — ^ be the eigenvalues of the symmetric tensor D 

D = AieiBj +\>e 2 ei + Aaeae 3 ™ 
Diffusion can be divided into three basic cases depending on the rank, of the representation tensor: 

1) Linear case diffusion is mainly in the direction corresponding to the largest eigenvalue. 
D ~ AiDf = A 1 e,e, J . 

(8) 

2) Planar case (Ai — \-> 3> diffusion is restricted to a plane spanned by the two eigenvectors 
corresponding to the two largest eigenvalues. 

D ~ 2AiD p = Ai(gjGj + e->e, ). 

(9) 

3) Spherical case ^ L — A^ — ^ ^isotropic diffusion: 
D ~ 3AiD a = Ai(Mf + e,kt + Mj}. 

....(10) 

In general, the diffusion tensor D will be a combination of these cases. Expanding the diffusion tensor using 
these cases as a basis gives: 

D = Aieiei + A^ei + A 3 e 3 e 3 

= (Ai k)h&i + {h A 3 )(eiej + e 2 ei) 
+ A 3 (M] +eieij + e 3 e] ) 

= (A, A,)D ( 4 (Aj A 3 )D p ■ A ;! D B 
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where f 1 - > 3 t 2 3 ■' and ' :! are the coordinates of D in the tensor basis ' Uf ■ D- D- .A similar 
tensor shape analysis has proven to be useful in a number of computer vision applications. As described,the 
relationships between the eigenvalues of the diffusion tensor can be used for classification of the diffusion 
tensor according to geometrically meaningful criteria. By using the coordinates of the tensor in our new basis 
measures are obtained of how close the diffusion tensor is to the generic cases of line,plane and sphere. The 
generic shape of a tensor is obtained by normalizing with a magnitude measure of the diffusion. Here we define 
this magnitude as the largest eigenvalues of the tensor.This gives for the linear, planar and spherical measures: 



A, A, 



An anisotropy measure describing the deviation from the spherical case is achieved as follows: 



= c t + c„ = 1 



...(13) 



7.4. Visualization of Diffusion Tensors 

A 3D diffusion tensor can be visualized using an ellipsoid where the principal axes correspond to the tensor's 
eigenvector system. However,it is difficult to distinguish between an edge -on, flat ellipsoid and an oblong one 
using the surface shading information. Similar ambiguity exists between a face-on,flat ellipsoid and a sphere. We 
propose two techniques for the visualization of tensor fields that overcome the problems with ellipsoids. We 
compares the ellipsoidal representation of a tensor with a composite shape whose linear,planar,and spherical 



components are scaled accordi 
moQt ,„ roo e i= r\,, and 



. andr 



can be used for visualiz 



(cl, cp, and cs).Additionally,coloring based on the shape 
i of shape. 



VIII. CONCLUSION 

In this paper,we sketched some of the fundamental concepts of medical image processing.lt is 
important to emphasize that none of these problem areas has been satisfactorily solved,and all of the algorithms 
we have described are open to considerable improvement. In particular,segmentation remains a rather ad hoc 
procedure with the best results being obtained via interactive programs with considerable input from the 
user.Nevertheless,progress has been made in the field of automatic analysis of medical images over the last few 
years thanks to improvements in hardware, acquisition methods, signal processing techniques, and of course 
mathematics. Curvature driven flows have proven to be an excellent tool for a number of image processing tasks 
and have definitely had a major impact on the technology base. Several algorithms based on partial differential 
equation methods have been incorporated into clinical software and are available in open software packages 
such as the National Library of Medicine Insight Segmentation and Registration Toolkit (ITK), and the 3D 
Sheer of the Harvard Medical School [90]. These projects are very important in disseminating both standard and 
new mathematical methods in medical imaging to the broader community. 
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(c) Merging of active contours, (d) Steady state. 

Fig. 5 Myocardium segmentation in MRI heart image with two merging expanding conformal active contours. 

The mathematical challenges in medical imaging are still considerable and the necessary techniques 
touch just about every major branch of mathematics. In summary, we can use all the help we can get!We have 
proposed measures classifying diffusion tensors into three generic cases based on a tensor basis expansion. When 
applied to white matter the linear index shows uniformity of tract direction within a voxel while the anisotropic 
index quantifies the deviation from spatial homogeneity. The non-orthogonal tensor basis chosen is intuitively 
appealing since it is based on three simple, yet descriptive, geometrically meaningful cases. We have described 
how tensor diffusion data can be processed without reverting to the use of only scalar measures of the tensor 
data. By staying in the tensor domain, cleaning up of the data can be done meaningfully with simple methods 
such as smoothing. We discuss addition of tensors geometrically and argue that adding tensors and vectors are 
different in that tensor summation gives more than the "mean" event due to more degrees of freedom.By using 
the geometric diffusion measures on locally averaged tensors local directionality consistency can be determined 
(e.g. existence of larger fiber tracts). We have proposed that this averaging approach can be used to derive a 
tensor field that can be used to describe macrostructural features in the tensor diffusion data. The linear measure 
cl derived from the averaged tensor field can for example be used for quantitative evaluation of fiber tract 
organization. We also have described how non-linear operations can be used to remap the eigenvalues of the 
diffusion tensors and given a sketch of how this can be used for tracking white matter tracts. 




Fig. 6. Axial brain images showing the three geometrical measures and diffusion tensor maps with three different 
smoothing parameters. 

Top: shows the geometrical measures and the tensor map derived from the original data. 

Middle: shows the same measures derived from data that has been averaged with 9x9x3 Gaussian kernel. 

Bottom: from data averaged with a 15x15x5 Gaussian kernel. The rightmost column shows the tensors. The blue 
headless arrows represent the in-plane components of cieiThe out-of-plane components of n. ei(cl A el) are 
shown in colors ranging from green through yellow to red, with red indicating the highest value for this 
component. 
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Figure 8.3. Left: Diffusion tensors, weighted with their linear measure cl r;i ., from an axial slice of a human 
brain. Right: Averaged diffusion tensors using a 5x5x3 Gaussian kernel weighted with their linear measure ci - 
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